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Abstract 



Consider the dynamics of a two-qubit entangled system in the decoherence environment, we 
investigate the stability of pairwise entanglement under decoherence. We hnd that for different 
decoherence models, there exist some special class of entangled states of which the pairwise entan- 
glement is the most stable. The lifetime of the entanglement in these states is larger than other 
states with the same initial entanglement. In addition, we also investigate the dynamics of pairwise 
entanglement in the ground state of spin models such as Heisenberg and XXY models. 
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I. INTRODUCTION 



Quantum entanglement is nonlocal and excess-classical correlation between separate par- 
ties, which is a most important character of quantum mechanics Lots of interest has 
focused on the nature of entanglement and the structure of entangled states P]. Besides, 
entanglement is also the most central and indispensable resource in quantum information 



processing such as quantum computation |2|], quantum teleportation [3|, quantum dense 
coding U and quantum key distribution (sj . In a word, quantum entanglement is not only 
importance in theory but also in practical applications. 

One the other hand, it is well known that decoherence is a vital factor that should not 
be neglected in quantum information processing. The coupling between any quantum sys- 
tem and its environment is inevitable. Thus the entanglement will evidently be reduced and 
even disappear because of this system-environment coupling, e.g. in a large scale quantum 
computer or during the course of entanglement distribution via noisy channels. The stability 
of entariglement depends on the initial entangled system - its entanglement structure and 
its size p, y, y, We may look on entanglement as a bond between different qubits, just 
like the chemical bond between different atoms. People studied the behavior of chemical 
bonds in different environment to understand how chemical bonds are formed. Thus inves- 
tigating the entanglement dynamics of different types of entangled states in the decoherence 
environment may help to gain some insight into the properties of the decoherence and the 
entanglement, which will provide useful hints for maintaining entanglement. And what kind 
of entanglement bond is the most stable under different decoherence models is an interesting 
problem. 

In this paper, we investigate the evolution of pairwise entanglement for two-qubit entan- 
gled states in the decoherence model which is described by general Pauli channels. We use 
the concurrence of Woottersjll|. which is related straightforwardly to the entanglement of 
formation (EOF), as the measure of entanglement for two-qubit entangled states. The most 
interesting problem is that given some general decoherence model, what kind of entangled 
states can maintain entanglement best. We find that with the same initial entanglement, the 
lifetime of entanglement in some specific class of entangled states is the longest. For a special 
decoherence model, that is depolarizing channels, all pure states together with some mixed 
states, which we call as Decoherence Path States (DPS) is the most entanglement-stable. 
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We present the analytic dynamics of two-qubit entanglement for these special entangled 
states. Furthermore, we also study the stability of the nearest neighbor entanglement in the 
ground state of some spin models such as Heisenberg and XXY model. Coincident ally, the 
conclusion is that in some noise models, the nearest neighbor entanglement in the ground 
state is also the most stable, though it is not maximized |12,]. 

The paper is organized as follows. In Sec. II we introduce the entanglement measure of 
two-qubit entanglement and the decoherence model, which can be viewed as a completely 
positive map. In Sec. Ill we investigate the dynamics of two-qubit entanglement under the 
influence of decoherence and try to find the special entangled states of which the pairwise 
entanglement is the most stable. In Sec. VI we examined the evolution of entanglement for 
some specific and maybe important mixed states, e.g. the ground states of spin models and 
the maximally entangled mixed states etc. In Sec. IV conclusions and discussions, together 
with some interesting open questions are presented. 



II. ENTANGLEMENT MEASURE AND DECOHERENCE MODEL 



There have been a number of measures for two-qubit entanglement, such as the entan- 



glement of formation 
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14| . negativity Q| and relative entropy of entanglement jl6] etc. 
In this paper, we adopt the well-established measure of entanglement concurrence as the 
measure of two-qubit entanglement. Consider a general two-qubit state, the density matrix 



is p. Then its time-reversed matrix is defined as 

'P = {Cly® ay)p*{Oy® Uy) (l) 

The concurrence of p is given by 

C = maxjo, v^- V^- V^- (2) 
where AjS are the eigenvalues of pp in decreasing order. The corresponding entanglement of 



i{C) = h[^-±^^^E^] (3) 



formation can be evaluated as 
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where h (x) = — xlog2 x — (1 — x) log2(l — x) is the Shannon's entropy function. 

However, it is not a very simple task to calculate the concurrence of a two-qubit system in 
an analytic way. Here we adopt a new method of calculating the entanglement of formation 



and thus the concurrence, which is based on Lorentz singular- value decomposition 17[. For 
an arbitrary 2x2 state p, there exists a 4 x 4 matrix with elements Rij = Tr{pai ® aj). In 
the real i?-picture, the density matrix p can be written as 

3 



P 



(4) 



ij=0 



where {(jj} are the Pauli matrices. 



Lemma 1: The 4x4 matrix R can be decomposed as i? = LiSL^, with Li, L2 finite proper 
orthochronous Lorentz transformations given by Li = T{A (g) A*)T'', L2 = T{B ® B*)T\ 
(l ■ ■ 1 \ 



where T = ^ 



1 1 

i —i 



. The normal form S is either of real diagonal form S 



(5) 



\l . . -1/ 

diag[sQ, Si, S2, S3] with Sq > Si > S2 > \s^\, or of the form 

f a ■ ■ b \ 

■ d ■ 

■ ■ d 

\c ■ ■ b + c — a J 

with a, 6, c, d real. And the Lorentz singular values of the second normal form are given by 
[so, si, S2, S3] = [^J{a - b){a - c), A/(a - b){a - c), d, -rf]. 

Lemma 2: Given a state p and the corresponding matrix R, the concurrence of p is C = 
max {0, (—So + si + S2 — S3)/2} depending on the Lorentz singular values of R. And Sq — Si — 
-32 + 53 = min Tr{LiRL^), where Li, L2 are proper orthochronous Lorentz transformations. 

Li,L2 

Based on the above two useful lemmas, we can see that the concurrence of a given 
density matrix p is only determined by the Lorentz singular values of the corresponding 
i?-matrix, which are the only invariants of a state under determinant 1 SLOCC operations 



17[. In the following, we adopt this method for calculating the residual entanglement of 
an initial entangled state in the decoherence environment. It is shown that the infiuence of 
the environment on the pairwise entanglement is refiected by changing the Lorentz singular 
values. Before proceeding to the details, we first introduce the decoherence model generated 
by Pauli operators. 



For a general decoherence model, it can be denoted as a completely positive map with 
an operator-sum representation. The effect of the general Pauli channels on a qubit p is 
described as follows 

3 

^{pi) = PoP + PiC^iPCTj (6) 

i=l 

3 

where pj > 0, = 1, and cxj are Pauli operators. This decoherence model includes some 
representative noise channels. When pi = P2= Ps it is just the depolarizing channel, which 
describes the decoherence process related to the couplings of quantum system to the thermal 
reservoir in the large temperature limit 0]. And when pi = P2 = the noise model is the 
dephasing channel, without energy exchange between the system and the environment, and 
only lose phase information. These kinds of decoherence models are common in several 
physical systems. In the rest of this paper, we will investigate how the decoherence will 
influence the pairwise entanglement in details. 

III. DYNAMICS OF PAIRWISE ENTANGLEMENT UNDER DECOHERENCE 
A. General Pauli Channels 

We fist assume that each qubit is independently coupled to the environment. The envi- 
ronment is characterized by the noisy channels generated by Pauli operators as shown in Eq 
(6). The initial state p, associated with this state R = LiSL^, is an entangled two-qubit 
state. Then after some time, p will be transformed to another state p' with much less en- 
tanglement, due to the action of the noisy channels, that is p' = £i£2(p)- According to Eq 
(6), it can be obtained easily that 



3 

p' =J2 i^^ ® Nj)p{Mi ® Ar,)t (7) 

i,j=0 



where Mj = ^Jpl^^i and = y/pjcrj- Now we can transformed this map into the -R-picture. 
Denote the i?- matrix associated with the state p' as R', then 

3 3 3 

R'=Y1 Lm.RL% = LmJR{J2 Ll,) (8) 

ij=0 i=0 j=0 

where L^, and L^. are Lorentz transformations given by Lm^ = T{Mi M*)T^ and L^^ = 
T{Nj (g) N*)Tl 



3 3 

For simplification, we can introduce Li = Yli^Mi and L2 = Yli^Ny Therefore, the 

i=0 i=0 

state evolution under decoherence is simply characterized by i? = LiRL^ in the real R- 

3 

picture. For the decoherence model we discussed here, Li = ^ T{Mi ® M*)T'^, L2 = 

2 = 

3 

T{Nj^N*)T'' . After simple calculation, it can be seen that Li = L2 = diag[l, Qi, Q2, Q3], 

3=0 

where Qi = PQ + P1-P2-P3, Q2 = P0-P1+P2-P3 and Q3 = Po-pi-p2+P3- It is obvious 
that the action of the noisy channels on the entanglement can be viewed as shrinking the 
Lorentz singular values by the above three coefficients. 

If the initial entangled states are set as pure states, then according to the Schmidt 

n 

decomposition theorem '2], an arbitrary two-qubit pure state \Q) can be expressed as 
\Q) = Ai|0'l') + A2|1'0'), where here Ai and A2 are non-negative real numbers satisfying 
A^ + A2 = 1. That is there always exist local unitary operations U and V , which satisfy 
\n) = {U (g) V)\no), where l^o) = Ai|01) + A2IIO). Here |0), |1) are the +1, -1 eigenstates 
of the Pauli matrix. In the i?-picture, |f2o)(^o| corresponds to the matrix: 

/ 1 • • A^-A2\ 
2A1A2 

■ 2A1A2 

. • -1 / 

And the local unitary operations on \Qq) correspond to left and right multiplication of 
Rq with orthogonal matrices, therefore in the i?-picture, an arbitrary pure state \Q) {Q\ 
corresponds to the matrix: 

R = LuRoLl (10) 

1 • 1 • 

where Lu = { ■ Oi ) = T{U ® U*)T^ and = ( ■ O2 ) = T{y ® V*)T\ with d and 

O2 are real 3x3 orthogonal matrices with determinant 1. 

In the following we will investigate how the pairwise entanglement changes under deco- 
herence in the real i?-picture and try to find what kind of entangled states, with the same 
initial entanglement, can maintain entanglement best. 

We starting by considering the initial entangled state is in the Schmidt decomposition 
form l^o). As discussed above, due to the coupling between the system and environment. 



Rn 



(9) 
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|Qo) is transformed into another mixed states p'^ — £i£2(|^^o) (^o|) ■ In the it!-picture, this 
action can be expressed as i?Q = LiRqL^, that is: 

/ 1 • • 

2AiA2(5| 



(Ai - A?)Q3 \ 



V (A? - A^)Q3 



-Ql ) 



Therefore the concurrence of can be obtained easily according to lemma 2. 



C = max{0, 



-] 



(11) 



(12) 



(13) 



where Co = 2A1A2 is the initial entanglement. Actually, this result can also be obtained by 
the conventional way of calculating the concurrence. Note that 

(A . . b\ 

■ D C ■ 

■ C E ■ 
\B . -A) 

where A = {\ - gi)/4, B = {Q\- QDW^jl, C = (Q? + g^)AiA2/2, D = (1 + gi)/4 + 
Q-ii}<\ — ^2)72 and E = (1 + (5i)/4 — Qz{}\ — A2)/2. For a state having a density matrix of 
the above form, the concurrence is given by C — max{0, Ci, C2}, where C\ = 2(|i?| — sjDE) 
and C2 = 2(C — A). We note C2 is always larger than Ci, thus C — max{0, 2(C — A)} 
which agrees with the above result. 

For an arbitrary pure state jJl), the corresponding it! matrix is shown in Eq. (10). In the 
decoherence environment, |n) is transformed to p'. In the real i?-picture, the evolution of 
the state is described by 

= L/iRLj^ — LiLijRqLi^Li^ (14) 

We then multiply Lj) and Ly to R from left and right respectively, and get another R- 
matrix R" = L'iRqL^, where L\ = LjjLiLu and L'2 = LyLl^Lv- This corresponds to 
local unitary operations on the state p' in the p-picturc, thus the concurrence of R" is 
identical to the one of R' . As we have discussed above, the concurrence of a given R 
matrix is only determined by the Lorentz singular values, and the action of the decoherence 
on the entanglement can be viewed as shrinking the Lorentz singular values. If we look 
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on {so, si, S2, S3) as the components of a vector, then the influence of decoherence is just 
shrinking this vector according to the coefficients Qi, Q2 and Q3. This can be reflected by the 
residual entanglement shown in Eq. (12). The action of the noisy channels is characterized 
by the three shrinking coefficients. In addition, we note that L^,Ly G 5*0(3). Therefore 
the effects of Lu and Ly on Li and L2 respectively is changing the shrinking directions. 
Then if we order the coefficients {rii, n2, ris} — {Qi, (52, Qa} such that < n| < n|. Based 
on the above analysis, it can easily be seen that the maximum residual entanglement is 

C;^=ma.{0. ^°'"- + 'f +"^-' } (15) 

where Co is the initial entanglement. The above maximum residual entanglement can be 
achieved by appropriate local unitary operations. The corresponding initial pure states are 

p = {\n) = {U(^ U)\no) I Lj;LrLu = diag[l, ±nl/^, ±nl/„ ±nl]} (16) 

This set of pure states present those special states that are the most entanglement-stable. 
The minimum residual entanglement and the corresponding initial pure states can be also 
derived in a similar way. Given a specific decoherence model, the relation between Qi, Q2, 
and Qs is known, then the pure states in p can be written explicitly. For example, in the 
dephasing channels model, the coefficients are Qi — Q2 — Po — Ps and Q3 — 1. Therefore 
both \Qo) and {ax ax)\^o) belong to the set p. The other states contained in p can also 
be obtained easily according Eq. (16). 

Although all the pure states with the same entanglement is equivalent under local unitary 
operations, it can be seen from the above results that under decoherence the behavior of 
different pure states are not all the same. This reflect the properties of the decoherence model 
and its influence on the entanglement. Our results suggest that if the decoherence model is 
flxed, there exists a special class of pure states, with the same initial entanglement, which are 
more favorable for maintaining entanglement. This gives some useful hints for maintaining 
and distributing entanglement. For example, if we want to distribute an entangled pure 
state between two separate parties through noisy channels, it will be helpful to apply some 
local unitary operations beforehand to transform the entangled pure state into the form of 
the states in the above set p. 

For an initial general mixed entangled state p with the initial entanglement Co, to drive 
an analytic evolution equation for its entanglement maybe intractable. However, we derive 
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a upper bound of the residual entanglement in the following. This upper bound is the 
corresponding residual entanglement C'^ for the states in the set p with the same initial 
entanglement Cq. It has be shown in there exists an optimal decomposition p = 
J2iPi\Vi) iv^il^ such that Cj = Cq for each \{pi), with J2iPi — 1- Then under decoherence p is 
transformed into another state p' = J2iPi^'i-^2{\fi){fi\)- According to the convexity of the 
concurrence, we know that C < J2iPiCi- Since < C'^, it is obvious that C < C'^. This 
suggests that the state in the set p is the most entanglement-stable of all the states with the 
same initial entanglement, no matter pure or general mixed states. It is well known that, an 
arbitrary pure state can be transformed into the states, with the same entanglement, of the 
set p by local unitary operations. Thus in the sense discussed in this paper, the pairwise 
entanglement in pure states are more favorable for maintaining entanglement compared to 
the generic mixed states. This conclusion is valid for any decoherence model which can be 
verified from the above discussions. 

Up till now, we have investigate the entanglement dynamics of two-qubit states assum- 
ing that each qubit is independently coupled to the environment. In the following, we 
consider the situation that only one qubit is under decoherence and the initial state is a 
pure entangled state. We are interested in whether the lifetime of entanglement also de- 
pends on the initial entanglement in this situation. As discussed above, the entanglement 
dynamics for pure states with the same initial entanglement are also dependent on their 
forms. To be comparable, we set the initial pure state in the Schmidt decomposition form 
l^o) = Ai|01) + A2|10). If only the first qubit is under decoherence, then the two-qubit system 
becomes p' = ei{\Qo) {Qq\). In the i?-picture, the matrix corresponding to p' is R' = LiRq, 
that is 

f 1 • ■ Ai-A?\ 

2A1A2Q1 

2A1A2Q2 

V(A?-A2)Q3 • • -Q3 / 

Therefore the concurrence of p' is 

C = 2max{0, ^(|Qi - Qal - Qs - 1), ydQi + Q2I + Qa - 1)} (18) 

This result shows that the lifetime of the entanglement is independent on the initial entan- 
glement. No matter how much the initial two-qubit state is entangled, it becomes separable 



R' 



(17) 
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in a constant time. This somewhat interesting phenomena reflect that entanglement is some 
kind of nonlocal property. 



B. Depolarizing Channels 

In the above, we discuss the dynamics of entanglement under the noisy channels generated 
by Pauli matrices. When the parameters satisfy Pi = P2 = = p/4 and po = 1 — 3p/4, 
the decoherence model in Eq. (6) is the depolarizing channels. The depolarizing channels 
describe the system-environment coupling in the large temperature limit T ^ 00. It can 
be realized by random Von Neumann measurements. Again we assume that each qubit 
is independently coupled to the environment. The shrinking coefficients are Qi = Q2 = 
Qs = 1 — p. Taking into account of the strength of the system-environment coupling and 
the interaction time, we can write 1 — p{t) = e~^^ jTj]. Thus the corresponding matrices in 
the i?-picture are Li = L2 = diag[l,e~'^^,e~'^^,e~^^]. Given an arbitrary initial entangled 
pure state the corresponding Schmidt decomposition normal form is l^o) associated 
with the i?-matrix Rq. Then after time t, the residual entanglement is dependent on the 
R" = L1R0L2 . Note that LjjLiLu = Li and LyL2Lv = L2 here. Therefore the concurrence 
at time t is 



From the above evolution function of the pairwise entanglement, we can find that the 
residual entanglement at time t, only depends on its initial entanglement Cq. Thus for 
all two-qubit pure states coupled with the same depolarizing environment, the stability of 
the entanglement is only governed by their initial entanglement, although these pure states 
could be in different forms. In other words, all pure states are the most entanglement- 
stable, need not to apply local unitary operations beforehand. Recalling the above analysis 
in the -R-picture, the reason for this interesting result is that the shrinking of the Lorentz 
singular values of the associated i?-matrix, which is introduced by the depolarizing channels, 
is isotropic. Besides, it is also obviously that even for two generic mixed state, if there are 
LU equivalent then the dynamics of two-qubit entanglement are also equal. Furthermore, 
there are some special mixed entangled states which have the same entanglement dynamics 
as pure states, that is also the most entanglement-stable. We will discuss the situation of 
mixed states in the next section. In Fig. 1 we present a visual example by plotting the 




(19) 
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dynamics of residual entanglement for pure states and some other generic mixed states with 
the same initial entanglement. 




FIG. 1: (Color online) The dynamics of residual entanglement in a two-qubit system. The initial 
states are chosen as pure states and some other generic mixed states with the same initial entangle- 
ment. We set Co = 2/3. Residual entanglement C as function of time t. Pure state (Solid Curve); 
some other generic mixed states (Dashed and Dotted Curve). 

IV. EXAMPLES: MIXED STATES 
A. Decoherence Path States 

Definition: Given a decoherence model characterized by a completely positive map A, 
Decoherence Path States (DPS) are those transient states p obtained from the pure states, 
that is 3 \ip) which satisfies p = A{\ip){ip\). 

We consider the depolarizing channels, and the initial entangled state p is a decoherence 
path state with the initial entanglement Cq > 0. Thus there exists with the entanglement 
C and some time to that satisfy p = A(to){\ilj){ip\) and Cq = {C + |)e"^''*" — |. Then after 
some time t, the decoherence path state p evolves to another state p'{t) with the entanglement 
C'{t) = {C + i)e-2''(*o+*) - i. This can be simplified to C'{t) = max{0, (Cq + i)e-2«* - i}. 
It is obvious that this is the same with the dynamics of entanglement for pure states with 
the same initial entanglement Cq, as shown in Eq. (19). In fact, the familiar Werner states 
belong to the decoherence path states. 
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Therefore in the depolarizing channels not only all pure states but also some special mixed 
states have the same entanglement dynamics. In other words, the pairwise entanglement of 
the decoherence path states are also the most stable. We plot the dynamics of entanglement 
in decoherence path states with different initial entanglement, as depicted in Fig. 2. 

From Fig. 2 (a) it can be seen that the entanglement in the two-qubit system decreases 
with time due to its interaction with the decoherence environment. There exists some time 
Tc when its entanglement C {t) = for t ^ Tc. Thus Tc is the critical time when the 
system entanglement disappears, i.e. the two-qubit system becomes separable. We can 
easily calculate the critical time: 

T. ^ (20) 

The relation between and the initial entanglement Co is depicted in Fig. 2 (b). Certainly, 
the lifetime of entanglement is longer if the initial entanglement is larger. For the singlet 
state, the lifetime of entanglement is kTc = 0.549. 

(a) X (b) 




FIG. 2: (Color online) The dynamics of residual entanglement in a two-qubit system with initial 
decoherence path states, (a) Residual entanglement C as function of the initial entanglement Co 
and time t. (b) Critical time kTc as function of the initial entanglement Cq. 

However, we can not exclude the possibility that there exist non-DPS states with the same 
entanglement dynamics as DPS states. In the following, we can see a concrete example. But 
as a special class of mixed states, the decoherence path states of a given decoherence model 
are expected to exhibit some other interesting properties. 

B. Ground State of Spin Models 

In this section we will investigate the stability of the nearest neighbor entanglement of 
the ground states of spin models. The Hamiltonian of the translationally invariant XXZ 
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spin chain with periodic boundary condition is 



^ = + + ^^^^^i] (21) 

when 7 = 1 the above Hamiltonian represents the Heisenberg antiferromagnetic modeL The 



Pi,i+l 



(22) 



2l| are 



ground state Pj|l2l.ll8l| is translationally invariant and the total z component of spin is zero. 
Thus the reduced density matrix of the i and i + 1 site is: 

f u \ 

z* y 
\ V / 

When u = V and x = y ii corresponds to the ground state of the Heisenber g antife rromag- 
netic model. The density matrices of the maximally entangled mixed states 
also in the above form. In the general Pauli channels of (6), the state Pi/i+i evolve to 

(a e\ 

5 F 
F* C 
\E* D / 

where A = urjf + vrj^ + {x + y)rjiri2, B = xrjf + yrjl + {u + v)riiri2, C = yrjf + xr]2 + {u + v)riir]2, 
D = vrjl + urjl + (x + y)riir]2, E = {z + z*)r]sr]4 and F = zrjl + z*rjl with r]i = po + Ps, 
V2 = P1+P2, V3 = Po — P3 and 174 = pi—p2- We note that r/i > jr/sl and 772 > \ri4^\. In addition, 
xy > because Pi,i+i is positive. Thus \E\ < {BCY^"^. Therefore the concurrence of p'j j+i 
is given by C = 2max{0, |F| - {ADf/"^}. 

For the special Heisenberg antiferromagnetic model and the depolarizing channels, the 
residual entanglement is C = max{0,?7|Co — 2rjirj2}, where Co = 2(|2;| — u). Note that 

1/2}, 



Pi,i+1 



(23) 



l+e-''*)/2, r]2 = (l-e-'^*)/2 and 7/3 = e"'^*. Thus C = max{0, (Co+l/2)e 



-2Kt 



which is the same as the pure states. Therefore in the depolarizing channels, the nearest 
neighbor entanglement of the ground state of the Heisenberg antiferromagnetic model is the 
most stable, though it is not maximized This result study the entanglement of the 

ground states of spin models from a new point of view. In addition, we can verify that 
Pi,i+i does not belong to the decoherence path states. This suggest that in this certain 
decoherence model, several mixed states other than decoherence path states are also the 
most entanglement-stable. 
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V. CONCLUSIONS AND DISCUSSIONS 



In conclusion, we have investigate the entanglement dynamics of a two-qubit system 
under a general decoherence model, that is Pauli channels. Given a decoherence model, 
we find the special class of pure states that are the most entanglement-stable and present 
the analytic entanglement dynamics of these states. Since any pure states with the same 
entanglement are LU equivalent, we show that pure states are more favorable for maintain- 
ing entanglement than general mixed states. Therefore in the situation of maintaining or 
distributing entanglement, it is helpful to using pure entangled states and to apply appro- 
priate local unitary operations beforehand to transform the entangled states to the most 
entanglement-stable form. Particularly, we investigate a certain decoherence model i.e. the 
depolarizing channels. In this case, a special class of mixed states that is decoherence path 
states are as entanglement-stable as pure states. The familiar Werner states are indeed DPS 
states. In addition, we investigate the entanglement dynamics of some specific class of mixed 
states, such as the ground states of XXZ, Heisenberg antiferro magnetic spin model and the 
maximally entangled mixed states. It is shown that in the depolarizing channels, the nearest 
neighbor entanglement of the ground state of Heisenberg antiferromagnetic spin model is 
coincident ally the most entanglement-stable. Another interesting result is that if only one 
qubit is coupled with the decoherence environment then the life time of entanglement is 
independent on the initial entanglement. This just refiects that entanglement is some kind 
of nonlocal property. 

For the decoherence model considered in this paper, we find the most entanglement-stable 
form of entangled states. The extension of this study to more general decoherence models is 
very meaningful, which is also related to the important work in j2^. Furthermore, we also 
introduce a special class mixed states i.e. decoherence path states and find that they are 
also the most entanglement-stable states in the depolarizing channels. However, they are 
expected to exhibit more interesting properties in general decoherence models. 
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